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FINDING  THE  INTERSECTION  OF  TWO  CONVEX  POLYHEDRA 

D.  E.  Muller*  and  F.  P.  Preparata^ 
Coordinated  Science  Laboratory 
University  of  Illinois  at  Urbana -Champaign 


Abstract 


Given  two  convex  polyhedra  In  three-dimensional  space,  we  develop 
an  algorithm  to  (1)  test  whether  their  Intersection  Is  empty,  and  (11) 

If  so  to  find  a separating  plane,  while  (111)  If  not  to  find  a point  In 
the  Intersection  and  explicitly  construct  their  Intersection  polyhedron. 
The  algorithm  runs  In  time  O(nlogn) , ;^ere  n Is  the  sum  of  the  numbers 
of  vertices  of  the  two  polyhedra.  The  part  of  the  algorithm  concerned 
with  (111)  (constructing  the  Intersection)  Is  based  upon  the  fact  that 
If  a point  In  the  Intersection  Is  known,  then  the  entire  Intersection  Is 
obtained  from  the  convex  hull  of  suitable  geometric  duals  of  the  two 
polyhedra  taken  with  respect  to  this  point.  


*Also,  Departments  of  Mathematics  and  of  Computer  Science. 

^Also,  Departments  of  Electrical  Engineering  and  of  Computer  Science 

This  work  was  supported  In  part  by  the  National  Science  Foundation  under 
Grant  MCS76-17321  and  In  part  by  the  Joint  Services  Electronics  Program 
under  Contract  DAAB-07-72-C-0259. 
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Finding  Ch«  InCnrMcClon  of  Cwo  convex  poljrhodre  In  Chroo-dlMnaloncl 

■p«eo  1*  « elaaslcal  problan  In  coapuCaClonal  gaonatry  [^3.  A alapla  but 

ClM>consuBlng  aolutlon  to  chla  problaa  la  known,  and  it  la  ao  trivial 

that  It  la  not  avan  worth  a lltaratura  rafaranca.  It  goaa  aa  follow#. 

Lat  <7  and  S ba  two  convex  poljrhadra,  teat  each  face  of  ^ agalnat  aach  face 

of  to  aaa  If  they  Intaraact;  If  no  Intaraactlon  la  found,  than  the 

Intaraactlon  of  the  two  polyhadra  la  a«pty,  otharwlaa  It  can  be  alaiply 

2 

eonatructad.  It  la  claar  that  auch  an  algorithm  could  uaa  0(n  ) oparatlona 
If  n la  tha  aum  of  tha  numbara  of  vartlcaa  of  ^ and  S. 

Tha  aaarch  for  a mora  afflclant  procadura  haa  In  tha  paat  mat  with  no 
auecaaa.  Savaral  facta,  howavar,  auggaatad  chat  a mora  afflclant  mathod 
ahould  axlat:  flrat.  It  la  wall>known  that  two  polygon#  in  tha  plana  can 
ba  Intaraactad  in  tlma  llnaar  In  tha  avn  of  chair  numbara  of  vartlcaa  [^3: 
aacond,  aavaral  analoglaa  axlat  batwaan  tha  plana  and  tha  apaca,  l.a., 
convex  hulla  of  n-polnt  aata  [33  and  maxima  of  seta  of  o vactora  [4]  can  ba 
found  in  time  O(nlogn)  In  both  two  and  three  dlmanalona.  In  aplCe  of  thla, 
no  genera Illation  of  tha  polygon  Intaraactlon  algorithm  haa  bean  found. 

In  apaclal  caaaa,  however,  batter  than  quadratic  tlma  mathoda  have  bean 
known  for  aoma  tlma.  Specifically,  after  tha  development  of  tha  Laa- 
Praparata  algorithm  for  locating  a point  in  a planar  aubdiviaion  Cs3,  it 

waa  raallcad^^^  that  tha  Intaraactlon  of  two  polyhadra  can  ba  found  In 
2 

tlma  O(nlog  n)  If  a vertex  of  one  polyhedron  llaa  Inalda  tha  other 
polyhedron. 

7T1 

Private  eovniunlcationa  batwaan  M.  I.  Shamoa  and  F.  P.  Praparata, 

May  19'6.  Sea  alao  C 1 3 p.  IbO. 
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Alternately  (aee  [1^1,  p.  162),  Shamos  conjectured  that  polyhedron 
interaectlon  could  be  obtained  as  a merge  step  of  a dlvlde-and-conquer 
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algorithm  for  the  intersection  of  half-spaces. 

In  this  paper  we  present  an  algorithm  for  solving  Che  problem  of 
intersecting  two  polyhedra  in  time  O(nlogn).  The  algorithm  tests  whether 
Che  intersection  is  empty  and,  if  not,  explicitly  constructs  it.  The 
approach  is  baaed  on  the  fact  that,  if  a point  p in  the  intersection  is 
known,  the  intersection  can  be  obtained  through  geometric  dualization. 
Specifically,  the  two  polyhedra  are  both  transformed  into  their  geometric 
duals  with  respect  to  p,  and  the  convex  hull  of  the  dualiwhich  can  be 
found  in  time  O(nlogn),  is  the  dual  of  the  intersection  polyhedron 
(Section  3).  When  such  a special  point  p is  not  available, by  deploying 
known  techniques  in  time  O(nlogn)  we  can  test  whether  the  inter- 
section is  nonempty  and,  if  so,  obtain  a point  in  the  intersection 
(Section  4).  The  algorithm  requires  that  each  polyhedron  be  represented 
by  a very  versatile  data  structure  called  the  doubly  connected  edge  list, 
which  can  be  obtained  from  a more  conventional  representation  in  time 
linear  in  the  number  of  vertices  (Section  2). 
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2.  Derivation  of  « Doubly  Connected  Edge  Llat  for  « Planar  Graph 

Let  V ■ {v. ] and  E • [e, ,...,e  } be  the  sets  of  vertices  and 
in  La 

edges  respectively,  of  a planar  graph  embedded  In  the  plane  without 
crossing  edges.  We  assume  that  (V,E)  be  represented  as  follows.  To 
vertex  Vj€  V there  corresponds  cell  rCj]  of  an  array  H[l:n3,  which 
contains  a pointer  to  the  first  term  of  the  cyclic  list  of  the  edges 
Incident  on  Vj,  arranged  In  the  order  in  which  they  appear  as  one  proceeds 
counterclockwise  around  v^.  The  latter  lists  are  realized  by  means  of  two 
arrays  VERTEX[l:na]  and  NEET[l:2a]  so  that  (VERTExCll ,NEXT[i])  Is  the  format 
of  the  list  nodes.  This  representation  of  the  graph  (V,E)  Is  precisely 
the  one  obtained  by  the  algorithm  of  Preparata  and  Hong  [ 3]  which  con- 
structs the  convex  hull  of  a set  of  points  In  three  dimensions:  Indeed  the 
surface  of  a convex  polyhedron  Is  topologically  a planar  graph.  We  shall 
call  this  collection  of  lists  the  vertex-to-edge  representation  of  a planar 
graph. 

Although  the  vertex-to-edge  list  is  one  of  the  most  commonly  used 
representations  for  a planar  graph.  It  has  the  disadvantage  that  the  dual 
graph,  l.e.,  the  eec  of  faces  of  the  original  graph.  Is  not  resdlly  avallsble. 
For  this  one  would  have  to  develop  the  face-to-edge  representation  of 
the  original  graph  In  which  each  face  refers  to  a cyclically  ordered  list 
of  edges  which  enclose  It. 

A more  convenient  representation  for  this  purpose  Is  one  which  we  shall 
call  the  doubly  connected  edge  list  (DCEL) , from  which  we  can  obtain 
equally  easily  Information  either  about  the  edges  incident  on  a vertex 
or  the  edges  enclosing  a face.  We  will  now  describe  the  DCEL  and  give  In  an 
appendix  the  algorithm  by  which  to  obtain  It  from  the  more  conventional 
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v«rC«x-to-«<lg«  rcpreacncacion  of  Che  graph  in  Clme  proporClonal  Co  Che 
nuaber  n of  verCleea. 

The  aaln  coaponenC  of  Che  OCEL  of  a planar  graph  (V,E)  la  Che  edge 
node.  There  la  a one-Co-one  correapondence  becveen  edgea  and  edge  nodea, 
l.e.,  each  edge  la  repreaenCed  only  once,  whereaa  In  Che  verCex- Co-edge 
llac  each  edge  appeared  Cwice.  An  edge  node  eonalaca  of  four  Informaclon 
flelda  VI,  V2,  FI,  and  F2  and  Cvo  polnCer  fielda  PI  and  P2:  Chare fore  Che 
correapondlng  daCa  acniccure  la  eaally  Implemenced  wlch  alx  arrays  wich 
Che  aane  namea,  each  conalaClng  of  m cells.  The  meanings  of  Chese 
fields  are  as  follows.  The  field  VI  conCalns  Che  name  of  Che  verCex  which 
la  Che  origin  of  Che  edge  , whereaa  V2  conCalns  Che  Cennlnus;  In  Chls 
manner,  Che  edge  receives  a convenClonal  orlenCaClon.  The  field  FI  and 
F2  conCaln  Che  names  of  Che  faces  which  lie  respecClvely  on  Che  lefc  and 
on  Che  rlghC  of  Che  edge  orienCed  from  VI  Co  V2.  The  polnCer  PI  (P2) 
polncs  Co  Che  edge  node  conCalnlng  Che  flrsc  edge  encounCered  afCer 
VI (V2)  when  one  proceeds  counCerclockwlse  around  VI  (V2) . Names  of 
faces  and  verClces  siay  be  Caken  as  Incegers.  As  an  example,  a fragmenC  of 
a graph  and  Che  corresponding  fragmenC  of  Che  DCEL  Is  shown  In  figure  1. 

IC  la  now  easy  Co  see  how  che  edges  IncldenC  on  a given  verCex  or 
Che  edges  enclosing  a given  face  can  be  obCalned  from  che  DCEL.  If  Che 
graph  has  n verClcea  and  f faces,  we  can  assume  we  have  cwo  arrays 
HV[l:n3  and  HF[l:f]  of  headers  of  Che  verCex  and  face  Uses;  Chese  arrays 
can  be  filled  by  a scan  of  arrays  VI  and  FI  In  clme  0(n).  The 
following  a era IghC forward  procedure,  VERTEZ(J  ),  obCalns  Che  sequence  of 
edgea  IncldenC  on  v,  as  a sequence  of  addresses  scored  In  an  array  A. 
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Figure  1.  Illustration  of  the  DCEL. 


VERTEX(J) 

1.  begin  aCi]  **  a^  •-  a *-  HvCj],  i *-  2 


2. 

If  VlCa]  ■ j then  a 

•-  PlCa]  else  a •-  P2[a] 

3. 

While  a a^  do 

A. 

begin  a[i]  •"  a 

5. 

If  VlCa]  - j 

then  a *■  PlCa]  else  a •-  P2Ca] 

6. 

i - 1+1 

and 


end 

Clearly  VERTEZ(j)  runs  in  time  proportional  to  the  number  of  edges  incident 
on  v^.  Analogously,  we  can  develop  a procedure,  FACE(j),  which  obtains  the 
sequence  of  edges  enclosing  fj,  by  replacing  HV  and  VI  with  HF  and  FI, 
respectively,  in  the  above  procedure  VERTIX(J).  Notice  that  the  procedure 


VERTEX  traces  the  edges  counterclockwise  about  a vertex  while  FACE  traces 
them  clockwise  about  a face. 

In  an  appendix  to  this  paper  we  shall  show  that  the  DCEL  of  a planar 
graph  can  be  obtained  froa  its  vertex-to-edge  list  in  time  linear  in  the 


number  of  vertices. 


3 . Finding  the  Intersection  of  Tvo  Polyhedra  when  a Point  in  the 
Intersection  is  Known. 

Assume  that  the  two  convex  polyhedra  d and  3 in  three-dimensional 

Euclidean  space  are  represented  by  DCZL's.  They  are  taken  as  having 

vertex  secs  and  and  face  sets  and  t q respectively.  With  each 

vertex  v c V-  'J  V-  we  associate  its  three  Cartesian  coordinates 
u 3 

x^(v),  x^Cv),  x^(v) . With  each  face  f € U we  associate  a half-space 
bounded  by  its  plane  and  containing  the  corresponding  polyhedron.  The 
half-space  of  f is  described  by  the  inequality 


n^(f)x^  + n2(f)X2  + a2(f)x2  -tdCf)  ^ 0 (1) 

where  n^(f),  n, (f),  n^Cf)  and  d(f)  are  four  parameters  characteristic  of  f, 
normalized  so  that  d(f)  is  either  1,  0,  or  -1.  Our  objective  is  to  obtain 
a DCEL  for  the  intersection  of  d and  3 along  with  the  coordinates  of  its 
vertices  and  the  parameters  for  its  faces. 

In  this  portion  of  the  analysis  it  is  assianed  that  a point  in  the 

intersection  of  d and  3 is  known:  Che  problem  of  finding  such  a point 

will  be  discussed  later.  Without  loss  of  generality,  the  point  in  the 

intersection  will  be  taken  as  the  origin,  since  if  it  is  not  the  origin  a 

simple  translation  of  coordinates  will  make  it  the  origin.  We  begin  by 

assuming  that  the  origin  is  actually  in  the  interior  of  each  of  the  two 

polyhedra  and  hence  in  the  interior  of  their  intersection.  In  this  case 

each  face  f € F_  U F-  represents  a half-space  of  the  form 
Cl  <3 

n^(f)xj_  + ^ other  words,  the  'constaric  d(f)  = 1 

for  all  faces  f. 

Now,  for  each  of  the  polyhedra  d or  3 there  is  a correspon«iag  polyhedron 

(D)  (D)  (D) 

d or  3 , respectively,  which  we  shall  call  its  dual.  The  dual  d oi  d 

is  obtained  by  reinterpreting  the  coefficients  n^^Cf),  n^(fj,  n^(f)  of  each 


face  f 6 as  che  coordiiiaCes  of  a corTes350iiding  vertex  v,  of 

Conversely,  the  coordinates  x^(v) , X2(v) , x^Cv)  of  each  vertex  of  are 

reinterpreted  as  che  coefficients  of  a corresponding  face  of  This 

transfonnation  nay  be  regarded  as  a conventional  doalization  about  the 

unit  sphere  with  center  at  the  origin,  where  points  at  distance  A from  Che 

origin  are  transformed  into  planes  at  distance  1/A  from  the  origin  and 

vice  versa  ([7],  p.233).  A similar  procedure  allows  us  to  form  Che  dual  of  8. 

We  note  that  this  doalization  procedure  is  only  possible  because  the 
origin  is  in  che  interior  of  the  polyhedron,  if  this  case,  the  dual- is  also 
a convex  polyhedron  containing  che  origin.  If  che  origin  is  not  in  che  interior, 
then  some  of  the  inequalities  (I)  would  require  d(f)  = 0 or  -1,  and  we  ha%'e 
not  defined  dual  points  for  such  half-spaces. 

■LeC  vertex  sets  of  and  8 respectively.  It  is 

easily  seen  that  the  convex  hull  of  the  union  of  ^ and  ^ is  che  dual 

of  the  intersection  of  ^ and  <5.  Hence,  in  order  to  find  the  inCerseccicn 

of  <7  and  8 we  may  simply  use  the  algorithm  of  Preparata  and  Hong  [l]  to  find 

Che  convex  hull  of  the  set  of  vertices  7^  U 3^  in  time  O(alogn)  , and 

7 “ 

upon  taking  che  dual  of  the  result  we  obtain  the  desired  polyhedron. 


Now  Lee  us  assume  that  the  given  point  in  the  intersection  of  Cl  and  8, 
i.e.,  the  origin,  is  not  in  the  interior  of  their  intersection.  Then 
certain  faces  f € U F^  have  d(f)  ■ 0.  In  fact,  these  are  exactly  the 
faces  which  pass  through  the  origin.  Let  F'  be  the  set  of  such  faces. 

To  each  face  f € F'  there  is  a corresponding  inequality  of  the  form 

nj^(f')x^  + ^ 0,  (2) 

obtained  from  (1)  by  replacing  d(f')  by  0.  A point  x in  the  interior  of 
a n 9 must  strictly  satisfy  all  inequalites  of  type  (1)  with  f € F L F , 
that  is,  none  can  be  an  equality.  Such  a point  x exists  if  and  only  if 
all  inequalities  of  type  (2)  with  f € F'  can  be  satisfied  strictly  by  some 
point.  To  determine  whether  there  is  such  a point  we  first  fix  - 1 and 
write  the  strict  form  of  (2)  as 

n^(f')x^  + n2(f')x2  > - n3(f’).  (2') 

Here,  by  normalizing  the  coefficients,  we  can  take  - n3(f')  as  either  1,  0, 
or  -1 . 

The  question  of  whether  or  not  (2')  can  be  satisfied  for  every 
f € F'  is  a two-dimensional  problem  of  the  type  we  are  solving  here  lor 
the  three-dimensional  case.  The  inequalities  of  (2')  collectively 
represent  a two-dimensional  convex  set  which  can  be  found  ([l],  p.  158>  in 
time  O(nlogn).  Actually,  a faster  computation  of  this  convex  set  is  possible 
the  inequalities  of  (2')  can  be  partitioned  into  two  sets,  depending  upon 
which  of  the  two  polyhedra  they  pertain  to;  each  such  set  corresponds  to  a 
polygon  in  the  plane  x^  - 1,  and  the  desired  convex  set  is  the  intersection 
of  the  two  polygons.  It  is  known  that  finding  these  polygons  and  their 
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intersection  runs  in  time  at  most  proportional  to  n Cl3. 

If  no  solution  is  found  in  the  above  case,  the  case  ■ -1  must  also  be 
tried.  This  problem  is  similar  to  the  previous  one  except  that  -n^(f')  is 
replaced  by  n^(f')  in  all  the  inequalities  (2*). 

Let  us  suppose  that  we  have  been  able  to  strictly  satisfy  all  inequalities 
of  the  type  (2)  with  x^  > 1 or  -1  using  the  above  method.  Clearly,  they  will 
remain  satisfied  if  the  vector  x is  multipled  by  a positive  scalar.  To 
strictly  satisfy  all  the  remaining  inequalities  in  (1)  whose  right  hand  sides 
are  all  1,  we  simply  choose  such  a scalar  which  makes  all  the  left  hand 
sides  less  than  1.  The  resulting  point  is  in  the  interior  ot  d r\  B. 

Tf  it  is  Impossible  to  strictly  satisfy  all  the  inequalities  of  (2), 
then  any  one  which  cannot  be  strictly  satisfied,  say 


n^(f")x^  + n^Cf'Ox^  + - 0 

represents  a plane  through  the  origin  which  contains  the  intersection  of 
d n B.  Thus,  the  intersection  of  d and  B may  be  found  entirely  within  this 
plane.  'Hiis  problem  Is  analogous  to  the  one  discussed  before  and  can  be  solved, 
as  we  saw,  in  time  0(n). 


S 
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4.  Finding  « Point  in  the  Intersection  of  Two  Polyhedra 

In  the  preceding  section  we  have  shown  that  the  intersection  of  two 
convex  polyhedra  can  be  obtained  when  a point  in  the  intersection  is  known. 
Thus,  if  the  intersection  is  nonempty,  all  that  is  needed  is  to  find  one 
such  point.  The  objective  of  this  section  is  the  implementation  of  this 
task. 

Given  a convex  polyhedron  <2^  > plane  is  called  a plane  of  support  of 
Ct  if  it  has  at  least  one  point  in  eonmon  with  Ct  and  all  interior  points  of 
Cl  lie  on  one  side  of  the  plane.  Hereafter  we  shall  only  consider  planes 
of  support  parallel  to  the  x^>axis  and  briefly  refer  to  them  as  vertical. 

The  intersection  of  £7  with  ltd  vertical  planes  of  support  is,  in  general, 
an  annular  region  R(£7)  of  the  surface  £7  which, in  the  absence  of  degeneracies, 
reduces  to  a cycle  of  edges.  The  projection  of  R(£7)  on  the  (x^,X2)  plane  is 
a convex  polygon  £7*  (figure  2),  which  is  the  convex  hull  of  the  projections 
of  the  points  of  £7  on  this  plane. 


Figure  2.  A convex  polyhedron  £7,  the  annular  region  R(£7)  and  the  projection 
polygon  £7*. 


a 
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Th«  ragloa  R((^  is  easily  obtained  from  the  DCEL  description  of 
<7  as  follows.  For  any  face  f^  of  the  nonnal  to  f^  is  the  vector 
(nj^(f^)  .OjCf^)  ,nj(f^))  . It  is  perpendicular  to  f^^  and  points  toward  the 
interior  of  Ct.  Given  any  edge  e of  (2^  let  f^  and  f^  be  its  adjacent 
faces.  Then  e € R(^  if  and  only  if 


•n3(fj)  ^ 0. 


(3) 


Therefore,  we  begin  by  scanning  the  edge  set  of  7 until  we  find  an  edge 
e which  belongs  to  R((^ , by  verifying  condition  (3  ).  At  this  point,  we 
select  one  of  the  two  vertices  of  e,  call  it  v.  Among  the  edges  incident 
on  V there  are  either  one  or  two  new  edges,  different  from  e,  which 
belong  to  R(^  and  can  be  easily  found  by  applying  the  VERTEX  procedure 
described  earlier,  to  the  DCEL.  Thus  we  can  advance  in  the  construction 
of  , which  will  be  completed  upon  re-encountering  the  initial  edge  e. 
Once  R((7)  has  been  computed,  is  trivially  obtained.  All  of  these 
operations  can  be  carried  out  in  time  proportional  to  the  number  |v^|  of 
the  vertices  of  7.  Thus  we  have  the  first  steps  of  the  algorithm, 
given  polyhedra  7 and  .3  with  |v^|  + |v^|  • n: 

if  if 

Step  I.  Find  7 and  d . (This  step  runs  in  time  O(n).’) 


Step  2.  Using  Shamos-Hoey ' s polygon  intersection  algorithm  [l],  find 
* ★ 

Che  Intersection  of  U and  5 . If  Che  intersection  is  empty,  halt,  for 
u n /3  is  also  empty.  Else  let  p*  be  a point  in  the  intersection  of  7 and 
(This  step  runs  in  time  0(n),  according  to  L ll . ^ 

if  if  if  if 

Under  the  projection  of  to  .7  , p - ;Xj^(p  >,X2(p  '>>  is  in 


I 


13 

general  Che  Image  of  a vertical  segment  of  which  reduces  to  a single 
point  in  some  cases.  In  any  case,  Che  preimage  of  p*  in  Cl  is  easily  found 
in  time  0(n)  as  follows.  For  each  face  f € ^e  determine  the  x^-coordinate 
of  the  point  on  the  corresponding  plane  which  projects  to  p*;  specifically, 
this  x^-coordinate  is 

a(f)  » - (n^(f)x^(p*)  + n2(f)x2(p*)  + dCfjj/n^Cf ) . 

Let  a' » min  a(f)  and  a"  = max  a'(f).  Then  a' 
n^Cf)  < 0 n^Cf)  > 0 

and  a",  with  a'  ^ a(",  are  Che  x^-coordinates  of  the  extremes  of  the 

segment  which  is  the  preimage  of  p*  in  C7;  we  similarly  define  8'  and  0", 
with  S'  2;  3",  for  the  analogous  segment  in  S.  If  the  two  segments  overlap, 
then  any  point  in  their  common  portion  also  belongs  to  the  nonempty 
intersection  of  Cl  and  S.  Otherwise  assume,  without  loss  of  generality, 
that  a"  > 9 ' . Then  we  define  the  near-sides  of  Cl  and  8 to  be  the  sets  of 
faces  [fjlf^  is  a face  of  Cl,  n^(fj)  < O}  and  [g^lg^  is  a face  of  '3,n2(g^)  > o} , 
respectively.  Clearly  both  near-sides  are  obtained  in  time  0(  | + l^jl)  = 0(n) 

by  traversing,  in  a straightforward  manner,  the  DCEL  descriptions  of  and  8. 

By  projecting  the  near-sides  of  Cl  and  8 on  Che  (x^^.Xj) -plane  we  obtain  two 
planar  straight-line  graphs  (PSLG's)  and  G^,  with  respective  vertex 
sets  and  V^.  Thus  we  have: 

Step  3.  If  the  pre-images  of  p*  in  d and  3 under  an  x^-pro jection 
have  a common  intersection,  then  halt,  for  any  point  in  this  intersection 
is  internal  to  d 3.  Otherwise  obtain  Gj  and  G^^.  (The  pre-images  of  p* 
are  found  in  constant  time;  and  are  found  in  time  0(n).)  j | 

i 

I 

I 

Let  3 be  the  closed  domain  contained  in  the  intersection  of  d*  and  ,3*.  ’ ! 


For  each  point  u c 3 it  is  convenient  to  define  Che  function  5<u),  called 


1 
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X3-dlscance,  as  follows.  If  a(u)  and  are  the  x^-coordinates  of  the 

points  on  the  faces  of  the  near-sides  of  G and  3,  respectively,  which 
both  project  to  u,  then 

5(u)  ■ Qf(u)  - 3(u)  . 


Let  us  now  enalyre  Che  function  5(u)  defined  on  Imagine  superimposing 
and  to  create  a new  vertex,  conveniently  called  a pseudo-vertex, 
at  Che  intersection  of  each  edge  of  G^  with  an  edge  of  G^.  Denoting  by 
V*  the  sec  of  pseudo-vertices  thus  obtained,  we  can  define  a new  PSLG  G* 
with  vertex  set  U U V*.  The  vertices  of  Vj  and  will  be  called 
true  vertices.  Thus  the  domain  is  subdivided  into  regions  by  G*.  Notice 
that  Inside  any  region  of  G-  the  function  a(v)  is  linear  in  the  (x, ,x  ) 
coordinates  at  v;  similarly,  for  the  function  0(v)  inside  any  region  of  G^. 
Thus  in  any  region  induced  by  G*  in  the  function  i5(v)«cjr(v)  - 0(v)  is  linear 
in  X3^(v)  and  • Moreover,  a(v)  is  convex -downward  and  3(v)  is  convex- 

upward;  it  follows  that  5(v)  is  a convex-downward  function.  We  conclude 
that  the  minimum  of  S occurs  at  a vertex  of  G*.  Notice  that  |v*|,  and 
hence  ^ ,could  be  0(n  );  in  fact,  it  is  not  hard 

to  construct  two  planar  graphs,  each  with  v vertices,  so  that,  when 

2 

superimposed,  (v  - 1)  intersections  of  edges  are  obtained. 


Since,  by  hypothesis,  «(p*)  • a"  and  3(p*)  “3',  we  conclude  that 

5(p*)  » a"  ~ 0'  > 0,  It  follows  that  the  intersection,  of  G and  3 is 

nonempty  if  and  only  if,  for  some  v € 5(v)  ^ 0.  Therefore,  either  we 

find  one  such  point,  or  show  that  min  S(v)  > 0. 

y € 3 

To  this  end,  we  begin  by  evaluating  S at  true  vertices  of  G*  in  3. 
This  is  easily  done  if,  for  each  a € 7^  we  determine  the  region  r(a)  of  G^ 
to  which  a belongs.  If  r(a)  is  not  the  infinite  region  of  the  plane  in 
the  subdivision  Induced  by  G^,  then  r(a)  corresponds  to  a unique  face  of 
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the  polyhedron  3,  and  5(a)  is  easily  computed.  Similarly,  we  can  compute 
5(b)  for  each  b € in  A The  determination  of  r(a)  has  been  called  the 
location  of  point  a in  the  planar  subdivlsioa  induced  by  G^C5]  and  could 
be  done  one  point  at  a time.  However,  a faster  algorithm  has  been  recently 
developed  [6  ] for  collectively  locating  all  the  points  of  a set.  According 
to  this  technique  the  members  of  can  all  be  located  in  regions  of  the  planar 
subdivision  induced  by  Gq  in  time  0((c  |v^l+ jv^l) log jv^|) , for  some  constant  c, 
and,  reciprocally,  the  members  of  can  all  be  located  in  regions  of  the 
planar  subdivision  induced  by  in  time  0(( [v^l+c jv^j) loglv^j) . Therefore 
in  total  time  0((  Iv^l+lv^j)  loglv^l . jv^j)  - 0(  ( lv^l+ jv -j)  log(  • Iv^j) ) - O(nlogn) 
all  the  true  vertices  can  be  located.  Once  a vertex,  say, of  has  been 

located  in  G^,  the  x^-coordinates  of  its  preimages  in  and  ,3  are  obtained 
in  constant  time.  This  is  summarized  as  follows: 

Step  4.  Locate  each  true  vertex  of  G-  in  the  planar  subdivision  induced 
by  G^  and  vice  versa.  (This  can  be  done  in  time  O(nlogn>  using  the 
algorithm  of  [6j.)  If  there  are  no  true  vertices  in  go  to  Step  7, 

Else  evaluate  5 at  each  true  vertex  of  G . (This  can  be  done  in 
additional  time  0(n).) 

Suppose  at  first  that  there  are  true  vertices  in  and  assume  that 

for  some  true  vertex  v (say,  v ? VL)  we  have  5 (v)  < 0.  The  vertical 

C( 

line  through  v Intercepts  the  near-side  of  iT  in  a vertex  a and  the  near- 
side of  3 in  a point  b,  and  obviously  o(v)  » X2(a)  < x^(b)  - 3(v).  Thus, 
we  have  a point  p*  for  which  ^(p*)  > (p*)  and  a point  v for  which 
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a(v)  <P(v).  Consider  now  the  plane,  parallel  to  the  x^-axis  and 
containing  the  points  p*  and  v.  The  intersection  of  this  plane  with 
the  two  polyhedra  S and  ? is  shown  in  figure  3 (where  the  points  p' 
and  p"  have  been  deftnedK  By  convexity,  Che  segments  ap'  and  bp" 
are  entirely  contained  in  JT  and  5 respectively,  and  so  their  point  of 
intersection  q belongs  to  the  intersection  of  S and  S.  The  coordinates 
of  q are  thus  obtained  by  straightforward  calculations. 

3 


b 


figure  3.  Finding  a point  in  the  intersection  when  5(p*>  > 0 and  '(v)  < 0. 

Assume  next  chat  5 (v)  > 0 for  all  true  vertices  v € Vj.  U V',  and 
let  V*  be  a true  vertex  such  chat  5(v*)  ■ min  [5(v)|v  £ Vi.  U 
We  cyclically  test  each  of  the  edges  of  o7  incident  upon  v*  to  determine 
whether  the  function  5 decreases  as  one  moves  along  the  edge  away  from  v*. 

If  it  fails  to  decrease  for  all  edges  incident  upon  v*,  then  v*  is  an 
absolute  minimum  of  the  function  5.  Since  5 (v*)  > 0,  the  polyhedra  iT 
and  3 do  not  intersect. 
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Step  5.  Obtain  v*.  If  v*  is  an  absolute  minimum  and  5 (v*)  > 0, 
u^j^t  for  CT  3 “ 3 ) if  i (v*)  ^ Q but  v*  is  not  an  absolute  minimum,  go 

to  Step  6;  if  5 (v*)  S 0,  then  there  is  a point  q ^ ^2  ? 3,  obtained  as 
the  intersection  of  the  diagonals  of  the  trapezoid  formed  by  the 
x^-projection  pre-images  of  p*  and  v*  in  :7  and  3.  (All  of  this  work  can 

be  done  in  time  0(n).) 

The  remaining  case  is  when  5 (v*)  > 0 but  v*  is  not  an  absolute  minimum. 
Then  5 decreases  as  one  moves  along  at  least  one  of  the  edges  - call  it  e - 
incident  upon  v*.  and  on  this  edge  we  locate  a pseudo-vertex  p.  One  such 
pseudo-vertex  must  exist,  for  otherwise  the  minimality  of  5(v*)  would  be 
contradicted.  Let  r(v*)  be  the  region  of  to  which  v*  belongs  (known 
from  Step  4);  to  locate  pseudo-vertex  p,  we  cyclically  test  each  edge  of 
r(v*)  in  turn  and  find  the  one  which  intersects  e.^^  Clearly  5(p)  < 5(v*). 


Step  6.  Locate  a pseudo-vertex  p adjacent  to  v*,  such  that  5(p)  < 6(v*). 
If  5(p)  S 0,  then,  since  5 (v*)  > 0,  there  is  a point  q € ^ D 5,  which  may 
be  found  as  in  Step  5;  otherwise  go  to  Step  7. 

We  must  how  consider  two  cases.  The  first  is  when  there  are  no  true 
vertices  in  (Step  4);  then  the  boundaries  of  a*  and  3*  must  intersect, 
so  the  point  p*  may  be  chosen  at  an  intersection  of  these  boundaries  and 
Is  therefore  a pseudo-vertex.  The  second  case  is  when  5(p)  > 0 (Step  6). 

Both  these  cases  are  treated  by  using  an  algorithm,  called  the  wandering. 
algorithm,  which  wanders  among  the  pseudo-vertices  of  G*  and  which  uses 
at  most  0(n)  time.  Thus  we  have: 


(^^If  V*  happens  to  belong  to  more  than  one  region  of  Gq,  then  **^863  of 
all  such  regions  may  have  to  be  tested  to  find  the  unique  one  which  Intersects 
In  any  case,  the  number  of  such  tests  is  0(n). 


e. 
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Step  7.  If  the  test  of  Step  4 fails  use  p*,  while  if  the  test 

of  Step  6 fails  use  p,  as  the  starting  point  of  the  wandering  algorithm 

(to  be  described  below),  either  to  find  a pseudo-vertex  p such  that 

5 ("p)  < 0,  to  which  the  method  of  Step  5 can  be  applied,  or  find  a 

pseudo-vertex  p such  that  5(p  ) ■ min  .5(v).  (We  shall  show  below  that 
m m , 

V t V 

the  wandering  algorithm  runs  in  time  0(n).) 

Before  describing  the  wandering  algorithm,  we  observe  that  the  starting 

point  of  it  is  a pseudo-vertex,  either  p or  p*,  which  has  a smaller  value  of  5 

chan  any  true  vertex.  If  we  imagine,  for  purposes  of  proof,  a contour  line  of 

' passing  through  p or  p*  we  enclose  a region  c which  contains  a 

pseudo-vertex  p having  minimum  5(p  ).  We  note  that  must  be  convex, 
m m 

Also  let  E'  and  be  the  sets  of  edges  of  and  G„  respectively  which 

intersect  S.  Since  no  true  vertices  lie  in  R,  each  edge  in  EL  U E' 

u S 

must  separate  R into  two  convex  regions.  No  two  edges  in  E^  can 
intersect  in  R,  nor  can  two  edges  in  Ei^.  Also,  the  function  ^ is  convex 
downward  as  one  travels  along  any  edge  of  E^  U E^  and  its  minimum  must  lie 
somewhere  in  R,  because  the  boundary  of  R is  a contour  line  for  5.  We 
shall  call  this  point  on  an  edge  e € E^  U E^  where  5 has 

a minimum  value,  a minimum  point  of  the  edge  e.  It  is  unique  except  in 
degenerate  cases.  The  value  of  5 at  this  point  will  be  called  the 
minimum  value  of  the  edge  e and  denoted  bv  min(e). 


' I 
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We  assume  momentarily  that  the  faces  of  both  polyhedra  are  triangles 

and, if  they  are  not, chat  we  shall  triangulate  both  polyhedra.  Notice  chat 

in  the  triangulated  polyhedra  Che  number  of  faces  remains  less  than 

twice  the  number  n of  vertices  and  Che  number  of  edges  remains  less  than 

three  times  n,  so  they  both  remain  0(n)  . Each  pseudo-vertex  p'  in  9.  is 

Che  intersection  of  two  edges  e'  € E'  and  e'  € E'  and  is  therefore  shared 

a b c 

by  four  regions  in  G*;  the  union  of  these  four  regions  is  referred  to  as 

the  crown  of  p*  and  is  the  locus  of  the  points  which  can  be  reached  from  p' 

without  crossing  anv  edge.  Notice  that  e'  is  shared  by  two  triangular 

a 

faces  of  G^,  whose  union  is  a quadrilateral  region;  a similar  remark  holds 

u 

for  e ' . Thus  Che  crown  is  the  intersection  of  these  two  quadrilateral 

D 

regions,  and  the  crown  boundary  contains  either  8,  or  10,  or  12  pseudo- 
vertices  (see  figure  ^ a,  b,  c,  respectively).  The  fact  that  the  number 
of  crown  vertices  is  bounded  is  a consequence  of  Che  hypothesis  that  the 
polyhedra  have  been  triangulated. 


(a)  (b)  (c) 

Figure  Illustration  of  the  possible  cases  for  the  crown  of  a pseudo- 

vertex. 

Given  any  pseudo-vertex  p'  as  the  intersection  of  e^  and  e^,  the  pairs 
of  triangles  bordering  these  two  edges  can  be  obtained  in  constant  time 
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from  the  doubly-connected  edge  lists  describing  and  ? respectively.  Once 
these  triangles  are  available,  the  pseudo-vertices  in  the  crown  can  also  be 
obtained  in  time  bounded  by  a constant,  and  so  can  their  values  of  5.  We 
now  give  the 

Advancing  step  of  the  wandering  algorithm;  A pointer  is  moved  from  the 
current  pseudo-vertex  p'  to  a pseudo-vertex  p"  which  attains  the  mininjum 
value  of  ; among  all  pseudo-vertices  in  the  crown  of  p'. 

Of  course,  the  step  is  voided  and  the  algorithm  terminates  if  p' 

attains  the  minimum  value  of  5 along  the  edges  e'  and  e'  intersecting 

^ D ] 

! 

in  p':  in  this  case  if  5 is  positive  the  two  polyhedra  do  not  intersect  1 

(case  (iii)  in  Section  3).  In  the  other  case  (5  decreases  either  along 

or  e^)  the  advancing  step  is  effected,  and  in  actual  practice  can  be  j 

carried  out  without  exploring  the  entire  crown  of  p',  but  simply  following 
a path  of  edges  along  which  6 decreases. 

An  additional  algorithm  simplification  is  that,  as  w®  shall  show, 

■! 

polyhedra  C7  and  3 need  not  be  triangulated  before  applying  the  wandering 
algorithm.  In  fact,  only  those  faces  of  and  G,  will  be  triangulated 
which  are  actually  traversed  by  the  wandering  algorithm.  Specifically,  let 
p',  the  Intersection  of  e^  and  e^^,  be  the  current  pseudo-vertex  (see 
Referring  for  simplicity  only  to  polyhedron  <7,  let  f^^  and  f^  be  the 
faces  of  G^  sharing  e\  In  the  doubly  connected  edge  list  of 


I 


I 

1 


Figure  5 ) . | 
two  I 
can 


Figure  5.  Partial  triangulation  ot  <2. 
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obtain  in  constant  time  the  edges  e"  and  e"'  which  follow  e'  in  the  edge-sequences 

of  f,  and  f„,  respectively.  If  f.  is  not  a triangle,  we  connect  the  non- 
1 Z L 

overlapping  extremes  of  e"  and  e',  and  we  do  likewise  for  f _ . The 
introduction  of  any  such  new  edge  in  the  doubly -connected  edge  list  requires 
the  modification  of  two  pointers  and  the  use  of  two  other  cells  for  con- 
struction of  the  appropriate  record.  All  this  can  also  be  done  in  constant 
time.  Ve  conjecture  that  this  insertion  is  not  really  necessary,  but  the 
present  proof  on  the  time  performance  of  the  algorithm  depends  upon  it. 

Since  the  wandering  algorithm  moves  from  p'  to  p"  only  if  5(p")  < 5(p'), 

it  is  obvious  that  the  algorithm  will  terminate  at  a point  p such  that 

m 

5(p  ) is  the  minimum  value  of  5 for  all  pseudo-vertices  in  ft.  Even  thou^ 
m 

2 

the  total  number  of  pseudo -vertices  in  ft  could  be  0(n  ) , we  shall  now  prove 

that  the  number  of  advancing  steps  is  at  most  0(n) . 

Recall  that  for  an  edge  e in  either  or  G^,  min(e)  denotes  the  minimum 

value  of  5 on  e.  Let  pseudo-vertex  p'  be  the  intersection  of  e'  € Ei  and 

a u 

®b  ^ ^3’  define  m(p')  » max(min(e '^)  ,min(e|^ ) ) . Clearly  5(p')  ^ ra'(p), 

Letiraa  1;  Let  p'  be  a pseudo-vertex  in  ft;  if  m(p'')  * 5(p  ).  then 
— m 

5(p')  - m(p')  = 5(p  ). 

m 

Proof:  Let  us  assume  the  contrary  and  obtain  a contradiction.  In  figure  6 

let  a*  and  b'  represent  minimum  points  on  e^  and  e^  respectively  which  are 

nearest  to  p'.  By  our  asstimption,  6 (a')  » S(^')  “ 5(p  ) and  hence  by 

m 

convexity,  every  point  along  the  line  segment  I between  a'  and  b'  also  has 
this  same  6 value.  Let  a^  be  the  pseudo-vertex  closest  to  a'  in  the  portion 
of  e'  between  a'  and  p'  (possibly,  a.  and  p'  coincide).  The  line  segment  I 
crosses  a region  of  G*  bordering  with  a'a..  Since  the  value  of  5 is  linear 


within  this  region  and  it  achieves  the  minimum  value  5(p  ) at  an  interior 

m 
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Figure  6.  Illustration  for  the  proof  of  3i(p')  ” ^(Pq)  * ^(P')  “ “(p'). 

point,  it  must  have  this  value  throughout  the  entire  region.  Hence, 

5 (a,)  ■ 5(p  ),  contradicting  our  assumption  that  a'  is  the  nearest  minimum 
1 m 

point  to  p'  on  e'.  This  proves  6(p')  ■ 6(p  ) .Q 
a m 

Assuming  now  chat  5(p')  > 5(p^),  we  see  by  Lemma  1 chat  m(p')  > SCp^^^) 
the  wandering  algorithm  is  applied  at  p',  it  steps  to  a new  pseudo-vertex  p". 
Lemma  2 : m(p")  < m(p'). 

Proof : We  distinguish  two  cases: 

(1)  p*  and  p**  do  not  belong  to  the  same  edge.  Let  p'  be  the  inter- 
section of  e'  and  e’  and  let  p”  be  the  intersection  of  e"  and  eV  (see 
a b a D 

figure  7(a)).  Let  p^  be  the  pseudo-vertex  in  ft  defined  above. 

We  claim  that  a straight  line  from  p to  p"  cannot  intersect  Che  interior 

n 

of  the  region  f of  G*  to  whose  boundary  p'  and  p"  belong,  except  at  p". 

In  fact,  if  it  did,  any  such  point  of  intersection  would,  by  convexity, 


When 


have  a value  of  6 as  low  as  6(p").  Since  p"  is  a minimum  point  of  f,  this 

would  imply  that  all  the  points  of  f have  the  same  value  of  6,  contradicting 

6(p')  > 6(p").  As  a consequence,  either  e^  or  e^'  separates  p^  from  p',so 

p belongs  to  the  shaded  regions  in  figure  7(a). 
m 

Assume,  without  loss  of  generality,  that  p',  and  hence  e^,  is  separated 

from  p by  e".  Then,  since  e'  does  not  cross  e"  in  the  straight  line 
m a a a 

between  p^  and  the  minimum  point  of  e^  intersects  e”  in  a point  a".  By 
convexity,  min(e')  ^ 5(a"),  with  equality  occurring  only  if  min(e')  » 5(p  ). 

ft  ft  & 

Assuming  eqtiality,  since  we  have  seen  that  m(p')  > 5(p  ) and  we  have 

n 

min(e')  ■ 5(a")  * min(e")  » 6(p  ),we  obtain  m(p')  >min(e")*  Assuming  instead 

ft  ft  n ft 

that  mln(e^)  > 5 (a"),  since  by  definition  m(p')  ^ min(ej^)  and  6 (a")  > mia(e^), 
we  also  obtain  m(p')  > min(e")< 

ft 

Two  subcases  must  be  considered.  First,  assume  p',  and  hence  e^,  is  also 
separated  from  p^  by  e|^’.  Then  by  an  identical  argument  m(p')  > min(eJJ),  so 
m(p')  > m(p")  « max(min(e^)  ,mln(e;^')) . Second,  assisne  it  is  not,  as  shown  in 
figure  7(a).  We  now  show  that  mia(e^)  S min(e^)  thus  reaching  the  same  conclusion. 


In  face,  since  5 (p”)  is  the  (ntnimum  value  in  f,  the  minim  point  on 

e"  occurs  either  at  p"  or  along  e"  on  the  opposite  side  of  p"  from  f.  A 

straight  line  drawn  between  this  minimum  point  and  intersects  at  a 

point  b"  such  that  5 (b")  ^ min(e^) , by  the  convexity  argument  used  earlier. 

But  5(b")  i min(e||) , whence  min(e")  2 min(e^') , as  claimed. 

(2)  p*  and  p"  belong  to  the  same  edge.  Without  loss  of  generality,  let 

p’  be  the  intersection  of  and  and  let  p”  be  the  intersection  of  e^  and 

e^  (see  figure  7(b)).  By  the  convexity  arg\iment,  is  separated  from 

p by  e"  (i.e.,  p belongs  to  the  shaded  region).  As  in  case  (1),  we 
m a m 

can  show  that  m(p')  2:  min(e^)  > min(e”).  To  prove  that  min(e^)  < mtn(ep 

we  note  that  the  minimum  point  of  e"  must  be  p",  for  otherwise  p"  would  not 

a 

attain  the  minimum  of  5 in  the  crown  of  p*.  Thus 

m(p")  *=  max(min(ep,  min(e^))  » rain(ep  < m(p'),  as  claimed4^ 

Theorem;  The  number  of  advancing  steps  performed  by  the  wandering 
algorithm  is  0(n). 

Proof ; We  have  shown  that  as  the  wandering  algorithm  moves  from  one 
pseudo-vertex  p'  to  the  next,  the  value  of  m(p')  decreases  at  each 
step.  Each  value  of  m(p')  is  the  minimum  value  of  one  of  the  edges  in 
E^  U E^.  Hence,  the  number  of  distinct  values  which  m(p')  can  assume 

«-4  tj 

is  no  greater  than  |e^|  + which  is  0(n).  The  number  of  steps 

taken  by  the  algorithm  therefore  is  0(n).Q 

Since  the  time  taken  by  the  wandering  algorithm  is  0(n),  the  time 
taken  by  the  entire  algorithm  remains  O(nlogn). 


5.  An  Application;  Fiadlng  a Separating  Plane 


The  preceding  method  can  be  used  to  solve  efficiently  the  important  problem 
of  linear  separability  in  three  dimensions,  l.e.,  testing  whether  two  finite  sets 
of  points  A and  B are  separable  by  means  of  a plane,  and,  if  so,  finding  one 
such  plane. 

Since  two  finite  sets  of  points  are  linearly  separable  if  and  only  if  their 
convex  hulls  do  not  intersect  [S],  we  begin  by  obtaining  the  respective  convex 
hulls  of  the  sets  A and  B by  means  of  the  Preparata-Hong  algorithm  [3].  Letting 
(a|  -t-  |b|  > n,  this  task,  which  is  completed  in  time  O(nlogn),  yields  two 
convex  polyhedra  Cl  and  B such  that  |v^|  + |v^|  "Cs.  n.  We  now  apply  to  Cl  and  B 
the  algorithm  described  in  Section  4:  any  time  the  algorithm  declares  that  Cl 
and  B do  not  intersect,  we  construct  a separating  plane. 

We  now  recall  that  Q and  B are  fotmd  to  be  disjoint  in  three  exclusive 
cases,  already  referred  to  iq  Section  4; 

(i)  after  projecting  R((7)  and  R(/S)  on  the  plane  (x^,X2),  the  polygons 
dif  and  B>f  are  disjoint; 

(il)  after  evaluating  6 at  all  true  vertices  of  G*  we  find  that 

6(v*)  <■  min  6(v)  > 0 and  v*  is  an  absolute  minimum; 

v€v^UV^ 

(ill)  after  applying  the  wandering  algorithm  we  find  that  S(Pjq)  0. 

In  case  (i)  it  is  sufficient  to  find  a straight  line  I separating 
Qtt  and  since  a plane  containing  I and  perpendicular  to  the  plane  (x^,X2) 
separates  Cl  and  B,  The  line  I can  be  found  in  time  0(n)  by  an  obvious 
modification  of  the  Preparata-Hong  algorithm  for  planar  convex  hulls  (C3],p.  90). 

Cases  (ii)  and  (ill)  can  be  handled  jointly  by  the  following  considerations. 
Rather  than  constructing  one  separating  plane,  we  construct  a locus  of  separating 
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plftOM  and  make  a salectlon  in  this  locus.  Let  u ba  the  point  at  which  the 

algorithm  terminates  with  the  result  that  (2  and  B do  not  intersect;  obviously, 

either  u"V*oru>p  . Also  let  u'  and  u"  be  the  pre-images  of  u (with 

m 

respect  to  x^-projection)  in  Ct  and  .respectively.  Assuae  at  first  that 

u • V*  and,  without  loss  of  generality,  let  u*  be  a vertex  (in  V^).  Consider  the 

cycle  F of  i^e  faces  sharing  u* ; for  each  f € F,  Imagine  applying  the  vector 

(nj^(f),n2(f),n2(f))  ■ n(f)  to  the  origin;  recall  that  n(f)  is  normal  to  f and 

pointing  toward  the  exterior  of  (2,  Then  the  set  of  directions  ^(f)|f  € f} 

defines  a convex  cone  such  that  any  direction  internal  to  it  is  normal  to 

a supporting  plane  of  2.  Notice  now  that,  when  u * p^,  point  u'  belongs  to  some 

edge  e^  of  Q and  degenerates  into  a plane  wedge  delimited  by  the  normals  to 

the  two  faces  of  (2  which  share  e_ . 

a 

For  u"  the  convex  cone  Zq  is  analogously  defined,  with  the  only  modification 
that  the  directions  of  the  vectors  n(f)  are  reversed.  The  cone  can  assume  the  follow- 
ing forms:  if  u ■ v*,  then  is  either  nondegenerate,  or  a plane  wedge,  or 
a half-line,  depending  upon  whether  u"  in  B is  either  a vertex, or  a point  in  an 
edge,  or  a point  in  a face,  respectively;  if  u ■ p^,  then  Zq  is  a plane  wedge. 

The  solution  to  our  problem  is  Notice,  however,  that  this 

intersection  consists  of  a single  ray  in  the  following  two  cases:  (1)  u ••  p , 
in  which  case  they  ray  is  the  conmon  normal  to  the  edges  which  contain  u'  and  u" 
in  <2  and  /9,  respectively;  (2)  u ■ v*  and  u"  is  a point  in  a face  of  B^  in  which 
case  the  ray  is  the  normal  to  this  face.  In  the  remaining  cases  (u  ■ v*,  u'  is 
a vertex,  and  u"  is  either  a vertex  or  a point  in  an  edge)  we  first  find  a 
plane  which  intersects  Z^  in  a bounded  polygon;  this  can  be  done  in  time 
0(n)  as  follows.  For  each  face  f € f,  let  point  t(f)  be  the  terminus  of 
the  vector  n(t).  Let  f^^  and  f2  be  two  consecutive  faces  in  the  cycle  F. 

We  consider  the  set  of  planes  determined  by  triples  of  points  (t(£j^) ,t(f2) ,t(fj)) 
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with  J ^ 1,2  and  € F.  For  each  such  plane,  ve  call  positive  the  half- 
plane bounded  by  the  straight  line  I through  t(£^)  and  tCf^)  and  containing 
t(£j).  The  half-planes  of  this  set  have  line  I in  common  and  are  comprised 
between  two  extreme  ones,  one  of  which  intersects  all  the  edges  of  C^: 
the  latter  defines  our  desired  plane.  Next  we  intersect  with  this  plane  and 
obtain  either  a polygon  or  a straight-line  segment:  in  any  case  the  problem  is 
reduced  to  finding  the  intersection  of  two  plane  polygons,  which  can  be  solved 
in  time  0(n)  [1].  This  enables  us  to  find  a vector  orthogonal  to  a 
separating  plane;  the  constniction  is  completed  by  requiring  that  the  plane 
contain  a point  internal  to  the  segment  u'u". 

Thus,  we  conclude  that  the  construction  of  a separating  plane  of  two 
three-dimensional  sets  of  points,  if  it  exists,  can  be  effected  in  time 


0<nlogn) 
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Appendix 

As  we  said  in  Section  2,  the  vertex- to-edge  list  of  a planar  graph 
is  a collection  of  edge  lists,  referred  to  as  input  edge  lists,  stored 
in  arrays  H[l:n],  VERTEX[l:2n],  and  NEXT[l:2n].  In  the  DCEL,  we  can 
identify  n cycles  of  edges  around  a vertex,  called  vertex  cycles,  and 
f cycles  of  edges  around  a face,  called  face  cycles.  The  construction  of 
the  DCEL  is  carried  out  in  two  phases.  In  the  first  phase,  we  fill  the 
arrays  VI,  V2,  PI,  and  P2,  hereby  constructing  the  vertex  cycles.  In  the 
second  phase  we  generate  the  names  of  the  faces  and  fill  the  arrays  Fl  and 
F2,  hereby  constructing  the  face  cycles. 

Informally,  phase-I  of  the  algorithm  works  as  follows.  The  input 

edge  1'*  ts  are  scanned  one  at  a time,  in  the  order  v ,v-,...,v  . While 

1 2 n 

scanning  the  input  edge  list  of  v^  an  edge  (Vj,v^)  is  entered  into  the  DCEL 
only  if  i > J:  in  this  manner  we  ensure  that  each  edge  is  entered  only  once. 

Thus  any  edge  (Vj,v^)  with  h < j is  already  present  in  the  DCEL,  since  it 
was  entered  while  scanning  the  input  edge  list  of  v^,  earlier  in  the 
execution  of  the  algorithm.  All  that  is  needed  now  is  therefore  the 
realization  of  the  appropriate  linking  of  such  (Vj,Vj^)  into  its  position 
in  the  vertex  cycle  of  v^ . To  effect  it  we  must  determine  the  location  of 
(Vj,v^)  in  the  DCEL.  This  can  be  done  as  follows, with  additional  storage 
0(n).  Suppose  that,  while  scanning  the  input  edge  list  of  v^^,  the  edge 
(v^,Vj)  is  to  be  entered  (obviously  h < J).  This  edge  is  linked  permanently 
into  the  vertex  cycle  of  v^^  and  temporarily  into  a list  of  edges  of  the 
form  (v^,Vj),  with  r < j.  The  members  of  the  latter  list  referred  to  as 
the  temporary  list  of  Vj,  are  linked  in  reverse  order  to  that  of  their 
occurrence  during  the  execution  of  the  algorithm.  Thus  this  list  can  be  managed 
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with  only  one  pointer  stored  In  an  array  IAST[l:n].  With  these  provisions, 
the  location  of  Is  easily  obtained:  In  fact,  prior  to  linking  the 

vertex-cycle  of  v^  we  scan  the  temporary  list  of  v^  starting  from  IAStCj] 
and  store  the  location  of  Into  cell  B[h]  of  an  auxiliary  array 

B[l:n].  Notice  that  the  latter  array  Is  only  scratch  memory  and  will  be 
used  repeatedly  for  each  v^.  Therefore  the  additional  storage  needed  consists 
of  the  arrays  LAST  and  B,  both  of  size  0(n},  and  of  program  variables 
a^,  *Q»  t,  r,  /• 

We  can  now  give  the  algorithm. 


CONSTRUCT  VERTEX  CYCI£S 


1.  begin  a «-  1 


for  J 1 step  1 until  n ^ IAStCj]  ♦“  A (Comment;  Initialize  LAST) 
for  j •“  1 step  1 until  n do 
begin  t - LAST[j] 

While  / A do 
begin  p •-  VICjC] 

bCp]  •- 1 

t - P2[i] 


Comment:  Loop  5-8  fetches  the  locations  of  all 
edges  with  r < j by  scanning  the  temporary 

list  of  Vj  and  stores  them  Into  the  array  B.  llils 
step  is  obviously  void  for  J ■ 1. 

t - hCj], 

r - VEFTEXCt] 
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^ r > J then 


besia 


Vl[a]-J,  V2Ca]  - r 
HVCj]  - a^  •-  a,  u »-  1 
P2[a]  - LAST[r] 
IAST[r]  - a 


a ♦-  a+1 


end  Comnent; Steps  10-15  initialize  the  vertex  cycle  for 
else  HV[j]  •-  a^  •-  sCr],  u •-  2 

Comment;  Steps  8-17  process  the  first  member  (v^,Vj) 
of  the  input  edge  list  of  v^.  If  this  edge  was  not 
previously  encountered  steps  11-16  are  executed; 
specifically,  the  edge  is  entered  in  step  12.  Variable 
a^  is  used  to  denote  the  location  of  the  last  member 
of  the  vertex  cycle  being  constructed. 

While  NEXTCt]  /fc  H[j]  ^ 
begin  t *-  NEXT[t] 

r - VEKTEXCt] 

^ r > j then 

begin  VlCa]  - j,  V2Ca]  *-  r 
P2Ca]  - LASlCr] 
lASlCr]  •-  PuCa^]  •-  a, 
aQ  - a,  u - 1 
a •-  a+1 


begin  PuCa^]  - B[r] 

Hq  - B[r],  u *-  2 


PuCa^]  - HVCj] 


Connent;  Steps  18-29  complete  the  construction  of  the  vertex 
cycle  forVj.  Specifically,  loop  18-28  successively  processes 
the  edges  incident  on  v^  and  either  enters  them  into  the  vertex 
cycle  (Steps  21-26)  or  simply  links  them  into  it  (Steps  27-28). 
Step  29  closes  the  vertex  cycle. 


To  evaluate  the  running  time  of  the  algorithm  Just  described,  notice 
that  each  edge  is  processed  exactly  twice:  once  to  be  entered  into  a vertex 
cycle  and  into  a temporary  list,  the  second  time  to  be  linked  appropriately. 
Both  these  operations  take  constant  time,  and  since  the  number  n of  edges  is 
0(a),  0(a)  time  is  used  to  fill  the  arrays  VI,  V2,  PI,  and  F2. 

To  complete  the  construction  of  the  DCEL  we  must  construct  the  face 
cycles.  The  next  algorithm,  CONSTRUCT  FACE  CYCIES,  starts  from  the  partial 
DCEL  which  is  produced  by  the  CONSTRUCT  VERTEX  CYCLES  procedure.  The  algorithm 
will  scan  the  DCEL,  using  an  integer  a as  a counter.  If  FlCa]  and  F2[a]  have 
already  been  filled, it  advances  to  the  subseqxient  edge;  otherwise  it  generates 
the  name  of  new  face  (using  a counter  s)  and  traces  the  ‘«dges  enclosing  it 
filling  the  appropriate  F-fields.  The  algorithm  terminates  when  2m 
filling  operations  have  been  performed:  an  integer  k is  used  to  control 


this  event 
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CONSTRUCT  FACE  CYCI£S 


1.  begin  for  j •-  1 step  1 until  tn  do  Fl[j3  •-  F2[j]  •-  A 


a •-  s *-  k •-  1 


While  k ^ 2in  do 


begin  If  Fl[a]  A and  F2[a]  A then  a - a+1 

else  begin  If  FlCn]  ■ A then  u 1 else  u •"  2 

Fu[a]  •-  s,  c •-  Vu[a],  HF[s]  - a^  •-  a , k •-  k+1 
a •-  PuCa] 

While  a ^ Sq  do 

begin  If  Vl[a]  ■ c then  u *-  2 else  u •-  1 
Pu[a]  •*  s,  c •-  Vu[a],  k •-  k+1 
a •-  Pu[a] 


s ♦“  s+1 


Since  In  the  latter  algorithm  each  field  Fl[a]  or  F2[a]  is  being 
processed  at  most  twice  (once  to  be  filled  in  steps  6 or  10,  and  possibly 
once  to  be  just  inspected  in  step  4),  the  running  time  is  0(n).  This  and 
the  analogous  result  for  the  vertex  cycle  algorithm  substantiate  our  claim 
that  the  DCEL  can  be  obtained  in  time  0(n)  from  the  original  vertex-to-edge 


list. 


